Atomistic Boron-Doped Graphene Field Effect Transistors: 
A Route towards Unipolar Characteristics 
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We report fully quantum simulations of realistic models of boron-doped graphene-based field effect 
transistors, including atomistic details based on DFT calculations. We show that the self- consistent 
solution of the three-dimensional (3D) Poisson and Schrodinger equations with a representation 
in terms of a tight-binding Hamiltonian manages to accurately reproduce the DFT results for an 
isolated boron-doped graphene nanoribbon. Using a 3D Poisson/Schrodinger solver within the Non- 
Equilibrium Green's Functions (NEGF) formalism, self-consistent calculations of the gate-screened 
scattering potentials induced by the boron impurities have been performed, allowing the theoretical 
exploration of the tunability of transistor characteristics. The boron-doped graphene transistors 
are found to approach unipolar behavior as the boron concentration is increased, and by tuning 
the density of chemical dopants the electron-hole transport asymmetry can be finely adjusted. 
Correspondingly, the onset of a mobility gap in the device is observed. Although the computed 
asymmetries are not sufficient to warrant proper device operation, our results represent an initial 
step in the direction of improved transfer characteristics and, in particular, the developed simulation 
strategy is a powerful new tool for modeling doped graphene nanostructures. 



The discovery of graphene has opened a promising al- 
ternative to silicon-based electronics. Reported charge 
mobilities in undoped graphene layers are actually or- 
ders of magnitude larger than those measured in silicon, 
but the unfortunate zero-gap semiconductor nature of 
this material severely limits the achievable lon/IoS ra- 
tio^ (ratio — infinite for an ideal switch — of the current 
flowing through the device in the conducting state to that 
flowing in the nonconducting state) for graphene transis- 
tors, making them not yet able to compete with main- 
stream silicon technologies 1^ Another adverse feature of 



graphene-based devices is their ambipolar electrical be- 
havior, which precludes the development of a comple- 
mentary logic architecture. 

A way to induce an energy gap in graphene is to 
reduce its lateral dimension. By using e-beam litho- 
graphic techniques and oxygen plasma etching processes, 
graphene nanoribbons (GNRs) down to a width of 10 
nanometers can be successfully fabricatediiSiii The re- 
sulting energy gap, however, is not sufficient for the 
achievement of a practically useful lon/IoS ratio in a 
transistor structure. Chemical, in-solution techniques 
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have been developed that allow one to obtain nanorib- 
bons with a smaller widtbi^ii^ but are hard to integrate 
into a large-scale device production technology. A mixed 
approach has been followed by Wang and Dai,^** who use 
lithography to define the nanoribbons, whose width is 
then shrunk by means of gas-phase chemical etching. Al- 
though quite interesting, this procedure could have some 
limitation in terms of the achievable device density. Fur- 
thermore, energy bandgaps are very unstable with re- 
gards to edge reconstruction and defects fiSi thus prevent- 
ing the fabrication of graphene transistors with perfor- 
mances comparable with their silicon counterparts. In- 
deed, almost ideal graphene nanoribbons have been ob- 
tained by means of careful chemical synthesis?^ however 
such an approach cannot be applied to the fabrication of 
complex circuits on insulating substrates in a straightfor- 
ward way. Strong chemical damage of graphene could in 
principle be used to tune the current flow, but with the 
drawback of a significant decay of the conductance that 
leads to poor performance in terms of current driveJ^^— 

Some of us have recently proposed a new class of de- 
vices based on the chemical doping of graphene by boron 
or nitrogeui^"— Using first-principles and mesoscopic 
quantum transport calculations (within the Landauer- 
Biittiker formalism) some indications of different behav- 
iors for electrons and holes were predicted in ribbons with 
lateral sizes as large as 10 nm, with mobility gaps in 
the order of 1 eVj^i^"^ Here the mobility gaps are typ- 
ically defined by the corresponding energy window for 
which the conductance is smaller than a typical value of 
G = lO-^Go (Go = 2eV/i) (see Ref,^ for details). Re- 
cent experimental findings confirmed the possibility to 
include boron and nitrogen impurities in graphene^"— 
therefore suggesting doped graphene transistors as a gen- 
uine alternative to ultranarrow GNRs with a more con- 
trolled tuning of the conductance features. However, 
fully self-consistent simulations of doped graphene tran- 
sistors are mandatory to ascertain the real impact of 
graphene doping. 

We explore the effect, at room temperature, of boron 
doping in graphene-based field effect transistors using 
the device geometry illustrated in Fig. [TJ Self-consistent 
NEGF calculations are performed, allowing for a proper 
treatment of the screened impurity potential and a more 
robust assessment of the doping effects on the electron- 
hole transport asymmetry and on the onset of mobility 
gaps. 



I. RESULTS AND DISCUSSION 

Parameters for our calculations have been obtained 
from a first-principle approach, as we detail in the fol- 
lowing. 

As an initial step, first-principle density functional 
theory (DFT) calculations were undertaken for a single 
substitutional boron impurity in two-dimensional (2D) 
graphene and in GNRSj^i^ by means of the SIESTA 



code,^^ within the local density approximation and us- 
ing a double-C basis set. Following the methodology in 
Refsi^i^S, the onsite energies of the Pz orbitals for the 
(bulk) 2D case were extracted. Their variation with the 
distance from the boron atom is shown in Fig. [2] (inset). 
The boron atom acts as an electron acceptor: the re- 
sulting impurity potential is repulsive for electrons and 
the screening length is found to be of the order of a few 
angstroms (the screening length^ is defined as the dis- 
tance from the charged impurity at which the surround- 
ing mobile charges, rearranging under the electrostatic 
action of the impurity, reduce its potential by a factor 
equal to Euler's number e). 

The transport properties of "quasi-metallic" and semi- 
conducting armchair GNRs (following Ref.--, we refer 
to an armchair GNR with N dimers contained in its unit 
cell as an "N-aGNR" ) with a single boron impurity at dif- 
ferent sites across the ribbon width are then computed, 
using the full ab initio Hamiltonian, in the Landauer- 
Biittiker framework.?^ Boron impurities actually induce 
a strong hole backscattering, with much weaker effect on 
the electron side. A similar result was already obtained 
in boron-doped carbon nanotubesr^i^ and attributed to 
strong backscattering at the resonance energies of the 
quasibound states localized around the boron impurity. 
Due to the lack of the rotational symmetry present in 
carbon nanotubes, in GNRs the energy of the conduc- 
tance dip on the hole side depends on the position of the 
impurity with respect to the ribbon edge (as illustrated 
in Fig. [land Fig.©. 

In Ref.'^'^ it has been shown that these ab initio re- 
sults can be fairly well reproduced with a simpler tight- 
binding model with a single orbital per carbon and boron 
atom. For each atom, the onsite energy is taken as the 
impurity potential extracted from the ab initio calcula- 
tion. The onsite energy of the boron atom is the only 
parameter that has to be adjusted in order to achieve the 
best match between tight-binding and ab initio transmis- 
sion profiles. The success of such a simple tight-binding 
model in describing boron impurities originates from the 
preservation of the sp^ hybridization of graphene in the 
presence of substitutional boron. 

By using this tight-binding model, transport studies 
in large armchair GNRs with a random distribution of 
boron atoms, with concentrations between 0.02% and 
0.2%, have been performed in Refi^. The addition of 
the transmission dips induced by all impurities leads to 
a mobility gap in the valence band, i.e. an energy inter- 
val where conductance is nearly suppressed by scatter- 
ing, while transport in the conduction band is preserved 
(the opposite behavior would be found in the case of ni- 
trogen doping). This occurs for semiconducting as well 
as for "quasi-metallic" GNRs, independently of ribbon 
width, as long as the potential can be assumed to be 
constant along the channel. The mobility gap increases 
with boron concentration and ribbon length. YLeiM- sug- 
gested a possible exploitation of this phenomenon to ob- 
tain field-effect transistors (FET) with wide GNRs (> 10 
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nm), within the reach of conventional lithography. The 
absence of a significant bandgap in these ribbons could 
be compensated for by the mobility gap. However, the 
above mentioned study has been performed for ribbons 
without any gate-induced charge in the channel, since 
the impurity potential considered was the one extracted 
directly from an ab initio calculation without taking into 
account the effects of screening when a gate is applied. 




FIG. 1. Sketch of the simulated double gate GNR field ef- 
fect transistor (FET): red dots represent the substitutional 
boron atoms. In the inset, the longitudinal cross-section of 
the transistor is shown (the scale is different along the x and 
y directions, to make the drawing clearer). 

In order to simulate the behavior of a realistic tran- 
sistor structure, we need to perform a self-consistent cal- 
culation to include the effects of an applied gate voltage 
and of a non-zero drain-source bias on the charge density. 
This can be done, at a reasonable level of computational 
complexity, with a tight-binding model set up with pa- 
rameters properly extracted from the ab initio results. 
However, differently from Ref.^'^, the impurity potential 
obtained from ab initio approaches cannot be directly 
used in a self-consistent simulation, since it would be 
overscreened, as a result of the Poisson/Schrodinger iter- 
ations. The approach chosen here is to mimic DFT self- 
consistency with the tight-binding model and a proper 
distribution of fixed charges; then the final impurity po- 
tential arises from the tight-binding Poisson/Schrodinger 
self-consistency. 

In particular, as we will detail in the Methods section, 
the ab initio results can be reproduced considering a fixed 
charge +e (where e is the elementary charge) in corre- 
spondence of each carbon atom, and using the following 
tight-binding representation; null onsite energies are con- 
sidered in correspondence of all atoms, while, following 
Ref.'^^, the hopping parameter is equal to tp — —2.7 eV 
for all the pairs of nearest neighbor atoms, with the ex- 
ception of the edge dimers, for which a value of 1.12 tp is 
used instead. 

To validate this approach, a self-consistent Pois- 
son/Schrodinger calculation (described in the Methods 



section) is performed for an isolated graphene layer with a 
single impurity, as in the ab initio calculation. As shown 
in Fig. [21 (inset), the obtained self-consistent potential 
around the impurity is very close to the ab initio result. 
Also the ab initio transmission spectra are accurately re- 
produced for all positions of the impurity, keeping the 
boron onsite energy equal to zero and with no further pa- 
rameter adjustment. Transmission spectra comparisons 
are shown in Fig. [2l 
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FIG. 2. Transmission through an isolated 32-aGNR (3.81 nm 
wide) with one boron impurity placed at two different dis- 
tances from the ribbon edge: 0.5a and 8o, where a = 2.46 A 
is the hexagonal lattice constant. The bold lines represent 
the results with our self-consistent tight-binding model. The 
thin dotted line represents the ab initio results. Inset: com- 
parison between the results for the impurity potential in 2D 
graphene obtained with self-consistent tight-binding and ab 
imtio approaches. 



Before performing full NEGF calculations, the effect 
of a gate with a nonzero applied voltage was analyzed on 
short ribbons with a single impurity. The ribbon is sur- 
rounded by a double gate with a dielectric thickness of 
1 nm, as shown in Fig. [T] The source-drain voltage and 
the Fermi level are set to 0. The charge and potential 
on the ribbon adjust self-consistently as a function of the 
gate potential Vg- First, the calculation was performed 
at Vg = with vacuum chosen as gate dielectric. The ob- 
tained potential is very close to the result discussed above 
for the neutral isolated ribbon. This gate configuration 
does not induce any significant additional screening. A 
more realistic model is obtained by considering a silicon 
dioxide dielectric between the gate plates and the ribbon, 
with relative dielectric permittivity Sr — 3.9. However, 
thin vacuum layers are kept just below and above the 
ribbon, with a thickness of 0.3 nm, of the order of the 
calculated graphene/Si02 distance.-^S. The impurity po- 
tentials and the corresponding transmission spectra are 
shown in Fig. [31 for Vg = and Vq = —0.5 V. Screen- 
ing is enhanced at a negative gate voltage, due to the 
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accumulation of holes in the ribbon. This has a strong 
effect on transmission: the conductance dip is less pro- 
nounced and is shifted in energy, which could compromise 
the presence of a clear mobility gap in a transistor con- 
figuration. It was verified that the effect of the gates is 
much larger than that of the thin air layers, in terms of 
the global electrostatics, and, in particular, of the screen- 
ing. This is the reason why in the runs for the evaluation 
of the transistor characteristics we will not introduce the 
air layers, that would make convergence of the Poisson- 
Schrodinger iterations much slower. 




Energy (eV) 

FIG. 3. Influence of the gate voltage on the transmission 
spectrum of a 32-aGNR with one boron impurity placed at a 
distance 0.5a from the edge. The calculation was performed 
with silicon dioxide as gate dielectric and a thin vacuum layer 
surrounding the graphene sheet. Inset: influence of the gate 
voltage on the impurity potential. 



To assess the performance of complete boron-doped 
graphene transistors biased in a typical operating 
point, we have used the open-source code NanoTCAD 
ViDES.— In this code, the 3D Poisson equation is solved 
self-consistently with the Schrodinger equation (with 
open-boundary conditions), within the Non-Equilibrium 
Green's Function (NEGF) formalism)^ Schottky con- 
tacts have been modeled at the GNR ends with the phe- 
nomenological approach described in Ref.— , including 
self-energies that mimic metallic contacts. 

The considered device structure is reported in the in- 
set of Fig. [H The device is a double-gate FET and the 
channel is a 32-aGNR (3.81 nm wide) embedded in Si02 
dielectric. This is consistent with typical nanofabrication 
techniques, in which the dielectric is deposited after the 
creation of the drain and source contacts. Top and bot- 
tom field oxides are 1 nm thick, and the channel length 
is equal to 20 nm. In Fig. 21 we show the transfer charac- 
teristics of the GNR-FET for a drain-source voltage Vds 
= 0.1 V, for different boron doping concentrations (de- 
fined as the ratio of the number of boron atoms to the 
total number of atoms in the channel). Dopant atoms 



are randomly placed along the channel. Results in Fig. |4] 
refer to a single distribution of dopants. 
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FIG. 4. Transfer characteristics for two different boron doping 
concentrations (0.3% and 0.6%), for a given random distribu- 
tion of the dopant atoms. The transfer characteristic for the 
undoped device is also reported. 

As expected, transfer characteristics for the undoped 
device are symmetric with respect to Vgs = Vds/ 2 — 
0.05 V. As soon as some boron doping is considered, the 
p-branch of the /ds-Vgs curve is suppressed, as pre- 
viously observed also in carbon nanotubes.^3,iiii The p- 
branch suppression increases as the doping concentration 
is increased. From these results, it appears that with 
the considered boron doping concentrations, while still 
not obtaining a sufficient lon/Ioff ratio for logic appli- 
cations^ with the considered 32-aGNR (which, however, 
having a "quasi-metallic" behavior in the absence of dop- 
ing, represents, in this regard, a worst-case scenario), a 
clear onset of unipolarity occurs. Increasing boron dop- 
ing, while further reducing the current in the p-branch 
(due to holes), also leads to a degradation of the device 
performance above threshold. 

In Fig. [5{a)-(b), we show the transfer characteristics 
for 23 different distributions of doping atoms along the 
channel, for the same doping concentrations as in Fig. 2] 
The thick solid lines correspond to the average transfer 
characteristics computed over the considered statistics of 
the doping distribution. As can be seen, different dop- 
ing distributions lead to a large dispersion of the transfer 
characteristics, which could be reduced if several nanorib- 
bon devices were operated in parallel (as can be necessary 
in some cases in order to increase the on current). 

Nevertheless, we remark that the electron-hole asym- 
metry is clearly developing with chemical doping, and, 
based on the results of Ref.-^S, it can be strengthened by 
increasing the gate length up to few hundreds of nanome- 
ters (especially for the case with 0.3%). Indeed, longer 
channel lengths could result in cumulated contributions 
of multiple scattering effects and localization, which will 
even more inhomogeneously impact on transport char- 
acteristics (as shown in equilibrium situations''^). Due 
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to computational limitations, we are not reaching such 
channel lengths here. Computational hurdles, in terms 
of difficulties in achieving self-consistence, also limit the 
doping level that can currently be simulated with our 
method to about 0.6%. Thus, at present we cannot han- 
dle the extremely high doping concentrations, of the or- 
der of 10%, that have been reached in a recent experi- 
mentSS. reporting unexpectedly large gap opening effects. 
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FIG. 5. Transfer characteristics for a) 0.3% and b) 0.6% boron 
doping concentration, considering a statistics of 23 randomly 
distributed dopant atoms. The thick solid lines correspond to 
the average transfer characteristic. For comparison, we report 
also the transfer characteristic for the undoped device. 



II. CONCLUSION 

In conclusion, by coupling a DFT-based tight-binding 
parameterization (validated via a multi-scale approach 
based on ah initio techniques) with self-consistent NEGF 
quantum-mechanical transport simulations, we have 
found the onset of an electron/hole transport asymme- 
try and mobility gaps in boron-doped graphene nanorib- 
bon transistors, confirming that chemical doping could 
serve as a way to shape the current- voltage characteris- 
tics of such devices. This is just a first step, because the 
achieved lon/Ioff ratio, although for a worst-case condi- 
tion, is still far too small for practical applications. There 
is however ample room for improvement and several dif- 
ferent approaches could be used to selectively incorporate 
chemical entities into graphene, thus engineering sepa- 
rate chemically modified areas from a single graphene 
layer. A doping strategy could therefore be elaborated 
to implement pn junctions, logic gates and complex ar- 
chitectures directly patterned on the same underlying 
graphene layer. Conventional state-of-the-art patterning 
methods, such as e-beam lithography, could be exploited 
to design superimposed material architectures (with ei- 
ther metallic or other semiconducting materials) onto the 
selectively doped and functionalized graphene-based sub- 
strates. Although our transistor simulations have been 



performed for graphene nanoribbons with a width be- 
low 5 nm (due to the otherwise prohibitive computa- 
tional cost), we wish to point out that the scaling anal- 
ysis detailed in Refs,-^"^- suggests that our main finding 
(electron-hole transport asymmetry) should be confirmed 
for ribbon widths above 10 nm, that is within the reach of 
conventional fabrication lithography techniques. Finally, 
longer ribbon lengths and low temperature are also fac- 
tors which are capable of widening the mobility gapsi^ 



III. METHODS 

The simulation approach is based on the observation 
that the DFT results can be recovered performing a self- 
consistent Poisson/Schrddinger calculation, within the 
NEGF formalism, on the system, which is described us- 
ing a tight-binding model and a proper distribution of 
fixed charges. 

In detail, we found that, following this procedure, the 
tuning of a single parameter of the tight-binding model is 
sufficient to properly reproduce the ah initio impurity po- 
tential and transmission spectra. The method is inspired 
by what actually occurs in the DFT calculation. The to- 
tal charge of a carbon nucleus and of its core electrons is 
-|-4e, while it is -|-3e for boron. In the isolated atoms, this 
charge is compensated by 4 valence electrons for carbon 
and 3 valence electrons for boron. In graphene, 3 elec- 
trons per atom (including, in the case of doped graphene, 
the boron atoms) tend to hybridize in a "sp^" fashion, 
while the remaining valence electron on each carbon lies 
on the "tt" bonding subband formed by the Pz orbitals. 
Also the Pz orbital on the boron atom is not empty, since 
it tends to attract an electron from the neighboring car- 
bon atoms, and this determines the shape of the impu- 
rity potential shown in Fig. [21 In the single pz orbital 
tight-binding approach, the "sp^" valence electrons are 
included in the global ionic charge of the atoms, which 
becomes -l-e for carbon and for boron. Therefore this 
is the net charge "localized" on each atom that we have 
to sum to the negative charge of the Pz atomic orbitals 
considered in the simulation. As far as the tight-binding 
parameters are concerned, the onsite energy for carbon 
atoms is 0, while the hopping energy between nearest 
neighbors is set to tp — —2.7 eV. The hopping energy of 
edge dimers is adjusted to \.12tp in order to reproduce 
the small bandgap obtained in DFT for "quasi-metallic" 
N-aGNRs (TV = 3p + 2, where p is an integer) 

Only the onsite energy of boron may be further ad- 
justed, performing a self-consistent Poisson/Schrodinger 
calculation for an isolated graphene ribbon with a sin- 
gle boron atom and comparing the obtained results with 
those deriving from the ah initio simulation of the same 
structure. 

For the Poisson equation the same boundary condi- 
tions are used in the ah initio and in the tight-binding cal- 
culations: Neumann boundary conditions (i.e. derivative 
of the potential energy, and thus electric field, constant 
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and, in particular, zero) in the longitudinal direction, at 
the ends of the nanoribbon, and periodic in the transverse 
direction, with a separation of 2 nm between the adjacent 
edges of nanoribbon replicas (indeed, with a large enough 
separation between the replicas this is equivalent to a 
Neumann boundary condition). The gate stack has been 
treated considering a constant permittivity, the one of sil- 
icon oxide, between the nanoribbon and each gate, where 
a Dirichlet boundary condition is assumed. At each iter- 
ation, the eigenvalues of the Hamiltonian are computed 
and the Fermi level is adjusted to obtain charge neutral- 
ity. The local charge density is then calculated and the 
Poisson equation is solved. The obtained electrostatic 
potential is subtracted from the onsite energies and the 
process is repeated until self-consistency is achieved. For 
all of the meshes tested in this work, it was possible to 
fit the ab initio data simply by adjusting the boron on- 
site energy. A typical choice was to spread the charge of 
each atom in a spatial region of half- width acci the inter- 
atomic distance, and to choose a mesh interval of 0.1 nm, 
small enough to achieve convergence of the Poisson so- 
lution. This charge repartition around each atom was 
defined using a smooth envelope function. In this case, 
the best fitting was obtained for a value of the boron on- 



site energy of about zero, thus substantially coincident 
with that for carbon. 
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